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Abstract 

Ambit stochastics is the name for the theory and appHcations of ambit fields and ambit pro- 
cesses and constitutes a new research area in stochastics for tempo-spatial phenomena. This paper 
gives an overview of the main findings in ambit stochastics up to date and establishes new results 
on general properties of ambit fields. Moreover, it develops the concept of tempo-spatial stochas- 
tic volatility/intermittency within ambit fields. Various types of volatility modulation ranging 
from stochastic scaling of the amplitude, to stochastic time change and extended subordination 
of random measures and to probability and Levy mixing of volatility/intensity parameters will be 
developed. Important examples for concrete model specifications within the class of ambit fields 
are given. 
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1 Introduction 



Tempo-spatial stochastic models describe objects which are influenced both by time and location. 
They naturally arise in various applications such as agricultural and environmental studies, ecology, 
meteorology, geophysics, turbulence, biology, global economies and financial markets. Despite the 
fact that the aforementioned areas of application are very different in nature, they pose some common 
challenging mathematical and statistical problems. While there is a very comprehensive literature on 
both time series modelling, see e.g. Brockwell & Davis ( 2002[ l, [Hamilton| ( [1994^ , and also on mod- 



elling purely spatial phenomena, see e.g. Cressie ( 1993 1, tempo-spatial stochastic modelling has only 



recently become one of the most challenging research frontiers in modern probability and statistics, 
see [Cressie & Wikle ( |2011| ), [Finkenstadt et al. ( 2007| ) for textbook treatments. Advanced and novel 
methods from statistics, probability, and stochastic analysis are called for to address the difficulties 
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in constructing and estimating flexible and, at the same time, parsimoniously parametrised stochas- 
tic tempo-spatial models. There are various challenging issues which need to be addressed when 
dealing with tempo-spatial data, starting from data collection, model building, model estimation and 
selection, and model validation up to prediction. This paper focuses on building a flexible, dynamic 
tempo-spatial modelling framework, in which we develop the novel concept of tempo-spatial stochas- 
tic volatility/intermittency which allows one to model random volatility clusters and fluctuations both 
in time and in space. Note that intermittency is an alternative name for stochastic volatility, used 
in particular in turbulence. The presence of stochastic volatility is an empirical fact in a variety of 
scientific fields (including the ones mentioned above), see e.g. jAmiri ( 2009[ ), [Huang et aL (2011 



Shephard & Andersen (2009 1. Despite its ubiquitousness and importance, however, this important 



quantity has so far been often overlooked in the tempo-spatial literature. Possibly this is due to the 
fact that stochastic volatility induces high mathematical complexity which is challenging both in terms 
of model building as well as for model estimation. 

The concept of stochastic volatility needs to be defined with respect to a base model, which we 
introduce in the following. While the models should be the best realistic description of the underlying 
random phenomena, they also have to be treatable for further use in the design of controls, or risk 
evaluation, or planning of engineering equipment in the areas of application in which the tempo- 
spatial phenomenon is considered. Stochastic models for such tempo-spatial systems are typically 
formulated in terms of evolution equations, and they rely on the use of random fields. We focus on 
such random fields which are defined in terms of certain types of stochastic integrals with respect to 
random measures which can be regarded as a unifying framework which encompasses many of the 
traditional modelling classes. We will work with random fields denoted by ^((x) G M with ^^(x) := 
y (x, t), where t G R denotes the temporal parameter and x G M'^ denotes the spatial parameter, where 
d G ]Nq. Typically, we have d G {1, 2, 3} representing, for instance, longitude, latitude and height. 
Note that by choosing continuous parameters (x, t), we can later allow for considerable flexibility 
in the discretisation of the model (including in particular the possibility of irregularly spaced data). 
This maximises the potential for wide appUcations of the model. We expect that the random variables 
Y'j(x) and It'(x') are correlated as soon as the points (x, t) and (x', t') are "proximate" according to 
a suitable measure of distance. This idea can be formalised in terms of a set At (x) C R'* x R such 
that for all (x',t') G ^i(x), the random variables lt(x) and ^('(x') are correlated. The set A((x) 
is sometimes referred to as causality cone and more recently as ambit set, see |B amdorff -Nielsen 



& Schmiegel (2004i, describing the sphere of influence of the random variable lt(x). A concrete 



example of an ambit set would, for instance, be given by a light cone or a sound cone. 
As the base model for a tempo-spatial object, we choose 

Yt{x) = [ h{x,t;ts)Lidtds), (1) 

where L is an infinitely divisible, independently scattered random measure, i.e. a Levy basis. Under 
suitable regularity conditions, our base model ([T]) can be linked to solutions of certain types of stochas- 



tic partial differential equations, which are often used for tempo-spatial modelling, see Barndorff- 
Nielsen, Benth & Veraart|(j201lll for details. The kernel function /i:R'^xRxR'^xR-;>R needs 



to satisfy some integrability conditions to ensure the existence of the integral, which we will study in 
detail in Section |2] Note that the covariance structure of the base model is fully determined by the 
choice of the kernel function h and the set At{x.), see Barndorff -Nielsen et al. (2010). In particular, 
by choosing a certain bounded set At(x), one can easily construct models which induce a covariance 
structure with bounded support; such models are typically sought after in applications, which feature a 
certain decorrelation time and distance. Under suitable regularity assumptions on h and on At(x), see 
Bamdorff-Nielsen, Benth & Veraart (2011 1, the random field defined in ([!]) will be made stationary 



in time and homogeneous in space. It should be noted that in any concrete application, one needs to 
account for components in addition to the base model, such as a potential drift, trend and seasonality, 
and observation error on the data level. 
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One of the key questions we try to answer in this paper is the following one: How can stochastic 
volatility/intermittency be introduced in our base model ([T])? We propose four complementary meth- 
ods for tempo-spatial volatility modulation. First, stochastic volatility can be introduced by stochas- 
tically changing the amplitude of the Levy basis L. This can be achieved by adding a stochastic 
integrand to the base model. This method has frequently been used in the purely temporal case to 
account for volatility clusters. In the tempo-spatial case, one needs to establish a suitable stochastic 
integration theory, which allows for stochastic integrands in the form of random fields. Moreover, 
suitable models for tempo-spatial stochastic volatility fields need to be developed. 

In the purely temporal, examples typically used to model stochastic volatility are e.g. constant 
elasticity of variance processes, in particular, square root diffusions, see Cox ( |1975| ), [Cox et al .|(|1985|), 
Omstein-Uhlenbeck (OU) processes, see Uhlenbeck & Omstein ( 1930 1 and more recently [BarndorfP] 
[Nielsen & Shep hard| ( |200T] ), and supOU processes, see |Bamdorff-Nielsen| ( [2000[ ), |Bamdorff -Nielsen] 
[& Stelzer, ( ,2011j ). Second, stochastic volatility can be introduced by extended subordination of the 
Levy basis. This concept can be viewed as an extension of the concept of stochastic time change as 
developed by Bochner ( 1949 1, see also Veraart & Winkel ( 2010| ) for further references, to a tempo- 
spatial framework. Last, volatility modulation can be achieved by randomising a volatility/intensity 
parameter of the Levy basis L. Here we will study both probability mixing and the new concept of 



Levy mixing, which has recently been developed by Bamdorff-Nielsen, Perez- Abreu & Thorbj0msen 



( |2012 1. We will show how these two mixing concepts can be used to account for stochastic volatil- 
ity/intermittency. 

Altogether, this paper contributes to the area of ambit stochastics, which is the name for the theory 
and application of ambit fields and ambit processes. Ambit stochastics is a new field of mathemati- 



cal stochastics that has its origin in the study of turbulence, see e.g. Barndorff-Nielsen & Schmiegel 



(20041, but is in fact of broad applicability in science, technology and finance, in relation to mod- 
elling of spatio-temporal dynamic processes. E.g. important applications of ambit stochastics include 
modelling turbulence in physics, see e.g. Bamdorff-Nielsen & Schmiegel (2004[ [ 2009[ ), [Hedevang 



(2012 1, modelling tumour growth in biology, see Bamdorff-Nielsen & Schmiegel (2 007] l, ponsdottir 



et al. (2008 ), and applications in financial mathematics, see Bamdorff-Nielsen et al.| ( j20T0| |, |Bamdorff- 



Nielsen, Benth & Veraart (2012 1, Veraart & Veraart (2012 ). 



The outline for the remainder of this article is as follows. Section |2] reviews the concept of Levy 
bases and integration with respect to Levy bases, where the focus is on the integration theories devel- 
oped by Rajput & Rosinski (19891 and Walsh| ( [r986 1. Integrals with respect to Levy bases are then 
used to establish the notion for our base model ([T]l and for general ambit fields and ambit processes in 
Section [3] In addition to reviewing the general framework of ambit fields, we establish new smooth- 
ness and semimartingale conditions for ambit fields. An important sub-class of ambit fields - the so- 
called trawl processes, which constitute a class of stationary infinitely divisible stochastic processes - 
are then presented in Section[4] Section[5]focuses on volatility modulation and establishes four com- 
plementary concepts which can be used to model tempo-spatial stochastic volatility/intermittency: 
Stochastic scaling of the amplitude through a stochastic integrand, time change and extended subordi- 
nation of a random measure, and probability and Levy mixing. Finally, Section[6]concludes and gives 
an outlook on future research. 



2 Background 

Ambit fields and ambit processes are constructed from so-called Levy bases. We will now review the 
definition and key properties of such Levy bases and then describe how stochastic integrals can be 
defined with respect to Levy bases. 
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2.1 Background on Levy bases 



Our review is based on the work by Pedersen ( 2003 | l, Rajput & Rosinski ( 1989 1, where detailed proofs 
can be found. 

Throughout the paper, we denote by {Q, T, P) a probability space. Also, let {S, S, Leb) denote a 
Lebesgue-Borel space where S denotes a Borel set in for a A; G N, e.g. often we choose S = U.^. 
Moreover S = B{S) is the Borel cj-algebra on S and Leb denotes the Lebesgue measure. In addition, 
we define 

Bb{S) = {A£S : Leb{A) < oo}, 
which is the subset of S that contains sets which have bounded Lebesgue measure. Note that since 



Leb is (T-finite, we can deduce that S = a{Bh{S)), see Peccati & Taqqu (2008 p. 399). Also, the 
set Bt,{S) is closed under finite union, relative complementation, and countable intersection and is 
therefore a 5-ring. 

2.1.1 Random measures 

Random measures play a key role in ambit stochastics, hence we start off by recalling the definition 
of a (full) random measure. 

Definition 1. 1. By a random measure M on {S,S) we mean a collection of M- valued random 
variables {M {A) : A G Bb{S)} such that for any sequence Ai, A2, ■ ■ ■ of disjoint elements of 
Bb{S) satisfying U'^^^Aj G Bb{S) we have M (u°^i^j) = J2T=i ^ (^i) 

2. By a. full random measure M on {S,S) we mean a random object whose realisations are 
measures on {S,S) a.s.. 

Note here that a realisation of a random measure is in general not an ordinary signed measure 
since it does not necessarily have finite variation. That is why we also introduced the term of a full 
random measure. Other articles or textbooks would sometimes call the quantity we have defined as a 



random measure as random noise to stress that it might not be a (signed) measure, see Samorodnitsky 
& Taqqu ( [1994[ p. 118) for a discussion of this aspect. 



In some applications, we work with stationary random measures, which are defined as follows. 

Definition 2. A (full) random measure on S is said to be stationary if for any s G 5 and any finite 
collection Ai, A2, . . . , An of elements (of B{S)) of Bb{S) the random vector 

(M (^1 + s) , M (^2 + s) , . . . , M {An + s)) 

has the same law as (M (^1) , M {A2) ,...,M {An)). 

The above definition ensures that a random measure is stationary in all components. One could 
also study stationarity in the individual components separately. 

2.1.2 Levy bases 

In this paper, we work with a special class of random measures, called Levy bases. Before we can 
define them, we define independently scattered random measures. 

Definitions. A random measure M on {S,S) is independently scattered if for any sequence Ai, A2, ■ . ■ 
of disjoint elements of Bb{S), the random variables M{Ai), M{A2), ... are independent. 

Recall the definition of infinite divisibility of a distribution. 
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Definition 4. The law /i of a random variable on M is infinitely divisible (ID) if for any n G IN, there 
exists a law on M such that /x = //*", where ^*"^ denotes the n-fold convolution of /i„ with itself. 

In the following, we are interested in ID random measures, which we define now. 

Definition 5. A random measure M on {S, S) is said to be infinitely divisible if for any finite col- 
lection Ai,A2,...,An of elements of Bb{S) the random vector (A/ (Ai) , M {A2) ,. ..,M (An)) is 
infinitely divisible. 

Let us study one relevant example. 

Example 1. Assume that M is an absolutely continuous fij.ll random measure on (S, S) with a density 
m and suppose that the stochastic process {m {x)}^^g is non-negative and infinitely divisible. Then 
M is an infinitely divisible full random measure. 

Now we can give the definition of a Levy basis. 

Definition 6. 1. A Levy basis L on {S,S) is an independently scattered, infinitely divisible ran- 
dom measure. 

2. A homogeneous Levy basis on [S, S) is a stationary, independently scattered, infinitely divisible 
random measure. 

2.1.3 Levy-Khintchine formula and Levy-Ito decomposition 

Since a Levy basis L is ID, it has a Levy-Khintchine representation. I.e. let C G IS and A G Bb{S), 
then 

C{CtL(^)}=log(E(exp(iCL(^))) 

= iCa*{A) - \eb*{A) + ^ (e<- - 1 - iCrf[_i,i](x)) n{dx,A), 



where, according to Rajput & Rosinski ( 1989 Proposition 2.1 (a)), a* is a signed measure on Bb{S), 
b* is a measure on Bb{S), and n(-, •) is the generalised Levy measure, i.e. n{dx, A) is a Levy measure 
on M for fixed A G Bb{S) and a measure on Bb{S) for fixed dx. 

Next, we define the control measure as introduced in Rajput & Rosinski (1989} Proposition 2.1 
(c). Definition 2.2). 

Definition 7. Let L be a Levy basis with Levy-Khintchine representation Q. Then, the measure c is 
defined by 



c{B) = \a*\{B) + b*{B)+ [ mm{l, x^)n{dx,B), B £ Sb, 

Jr 



(3) 



where | • | denotes the total variation. The extension of the measure c to a a-finite measure on (5, S) 
is called the control measure measure of L. 

Based on the control measure we can now characterise the generalised Levy measure n further, 
see|Rajput & Rosinski|(|1989[ Lemma 2.3, Proposition 2.4). First of all, we define the Radon-Nikodym 



derivatives of the three components of c, which are given by 

/ \ da* 1 / \ db* /, N n{dx,-) 

a(z) = — (z), 6(z) = — (z), z.(dx,z) = ^^-^(z). (4) 

Hence, we have in particular that n{dx,dz) = i'{dx,z)c{dz). Without loss of generality we can 
assume that u{dx, z) is a Levy measure for each fixed z and hence we do so in the following. 
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Definitions. We call {a,b,i'(dx, ■),€) = {a{z),b{z),v'(dx,z),c{dz))zes '^characteristic quadruplet 
(CQ) associated with a Levy basis L on (5, S) provided the following conditions hold: 

1 . Both a and b are functions on 5, where b is restricted to be non-negative. 

2. For fixed z, ^{dx, z) is a Levy measure on M, and, for fixed dx, it is a measurable function on 
S. 

3. The control measure c is a measure on {S, S) such that Jg a{z)c{dz) is a (possibly signed) 
measure on {S, S), b{z)c{dz) is a measure on {S, S) and i^(dx, z)c{dz) is a Levy measure 
on M for fixed B ^S. 

We have seen that every Levy basis on (5, S) determines a CQ of the form (a, v{dx, •), c) = 
(a(z), fo(z), v{dx, z), c((iz))z;g5. And, conversely, every CQ satisfying the conditions in DefinitionjS] 
determines, in law, a Levy basis on (5, S) . 

In a next step, we relate the notion of Levy bases and CQs to the concept of Poisson random 
measures and their compensators. 

Definition 9. A Levy basis on {S, S) is dispersive if its control measure c satisfies c({z}) = for all 
z G 5. 

For a dispersive Levy basis L on {S, S) with characteristic quadruplet (a, 5, c) there is a modi- 
fication L* with the same characteristic quadruplet which has following Levy-Ito decomposition: 

L*{A) = a*{A) + W{A)+ f y{N - n){dy, A) + [ yN{dy,A), (5) 

for A G Bb{S) and for a Gaussian basis W (with characteristic quadruplet (0, b, 0, c), i.e. Ty(74) ~ 
N{0, J^b{z)c{dz))), and a Poisson basis (independent of W) with compensator n{dy;A) 



E{N{dy;A)} where n{dx,dz) = i/{dx,z)c{dz), cf. Pedersen (2003 1 and Barndorff -Nielsen & 
|Stelzer| ( [20TT] Theorem 2.2) 

It is also possible to write ([5]) in infinitesimal form by 

L*{dz) = a*{dz) + W{dz)+ f xN{dx,dz)+[ x{N - n){dx,dz). (6) 

J{\x\>l} J{\cc\<l} 

This is particularly useful in the context of the Levy-Khintchine representation, which can then also 
be expressed in infinitesimal form by 

C{C t L{dz)} = log (E(exp«i(dz))) 

= iCa*{dz)--C^b*{dz)+ [ (e*^^' - 1 - <xlr i dz) 

■'^ ' (7) 

= (^iCa{z) - ^C^b{z) + j {e'^'' - 1 - iCa;I[_i_i] (x)) v{dx,z)^ c{dz) 

= C{CXL'{z)]c{dz), CGM, 

where L'{z) denotes the Levy seed of L at z. Note that L'{z) is defined as the infinitely divisible 
random variable having Levy-Khintchine representation 



C{C \ L'{z)} = iCa{z) - ^eb{z) + ^ (e<^ - 1 - iCxI[_i,i] (x)) iy{dx, z). 



(8) 



Remark 1. We can associate a Levy process with any Levy seed. In particular, let L'{z) denote the 
Levy seed of L at z. Then, {L[{z))t denotes the Levy process generated by L'{z), which is defined as 



the Levy process whose law is determined by Li(z) ^= L'{z) 
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Definition 10. Let L denote a Levy basis on {S, S) with CQ given by (a, b, v{dx, •), c). 
L If v{dr, z) does not depend on z, we call L factorisable. 

2. If L is factorisable and if c is proportional to the Lebesgue measure and a(z) and b{z) do not 
depend on z, then L is called homogeneous. In that case we write c{dz) = vLeb{dz) = vdz for 
a positive constant v > and where Leb{-) denotes the Lebesgue measure. 

In order to simplify the exposition, we will throughout this paper assume that in the case of 
a homogeneous Levy basis the constant v is set to 1, i.e. the measure c is given by the Lebesgue 
measure. 

2.1.4 Examples of Levy bases 

Let us study some examples of Levy bases L on (M'^, i3(M'^)) with CQ (a, 6, v{dx, •), c). 

Example 2 (Gaussian Levy basis). When u{dx, z) = 0, then L constitutes a Gaussian Levy basis 
with L{A) ~ N (^f^a{z)c{dz), f^b{z)c{dz)^, for A G Bi,(M.^). If, in addition, L is homogeneous, 
then L{A) ~ N {aLeb{A), bLeb{A)). 

Example 3 (Poisson Levy basis). When c{dz) = dz and a = b = and v{dx;z) = X{z)6i{dx), 
where 5i denotes the Dirac measure with point mass at 1 and A(z) > is the intensity function, then 
L constitutes a Poisson Levy basis. If, in addition, L is factorisable, i.e. A does not depend on z, then 
L{A) ~ Poisson {XLeb{A)), for all A e Sf,(M^). 

Example 4 (Gamma Levy basis). Suppose that c{dz) = dz, a = b = and the (generalised) 
Levy measure is of the form v{dx:,z) = x~^e~"^'^^^dx, where a{z) > 0. In that case, we call the 
corresponding Levy basis L a gamma Levy basis. If, in addition, L is factorisable, i.e. the function a 
does not depend on the parameter z, then L{A) has a gamma law for all A G Bh{^^)- 

Example 5 (Inverse Gaussian Levy basis). Suppose that c{dz) = dz, a = b = and the (generalised) 
Levy measure is of the form v{dx]z) = x~^/^e~2T (^)^dx, where 7(z) > 0. Then we call the 
corresponding Levy basis L an inverse Gaussian Levy basis. If, in addition, L is factorisable, i.e. the 
function 7 does not depend on the parameter z, then L{A) has an inverse Gaussian law for all 
A G ^b(M*^). 

Example 6 (Levy process). Ifk = l, i.e. L is a Levy basis on M, then L{[0,t]) = Lt, t > is a Levy 
process. 



2.2 Integration concepts with respect to a Levy basis 



In order to build relevant models based on Levy bases, we need a suitable integration theory. In 
the following, we will briefly review the integration theory developed by Rajput & Rosinski (1989l 
and also the one by |Walsh| ( 1986|), and we refer to ,Bamdorff -Nielsen, Benth & Veraart, ( ,2011 ) for 



a detailed overview on integration concepts with respect to Levy bases, see also Dalang & Quer- 



Sardanyons ( |2011[ ) for a related review and Basse-O'Connor et al. (2012) for details on integration 



with respect to multiparameter processes with stationary increments. 



2.2. 1 The integration concept by [Rajput & Rosinski] ( |1989 1 

According to Rajput & Rosinski ( 19891 p.460), integration of suitable deterministic functions with 
respect to Levy bases can be defined as follows. First define an integral for simple functions: 
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Definition 11. Let L be a Levy basis on (5,5). Define a simple function on S, i.e. let / := 

Yyj=i^j^Aj, where Aj e Bb{S), for j = 1, ... ,n, are disjoint. Then, one defines the integral, 
for every A G a{Bb{S)) = S, by 



fdL ■.= ^XjL{AnAj). 
7=1 



The integral for general measurable functions can be derived by a limit argument. 

Definition 12. Let L be a Levy basis on (5, S). A measurable function / : {S, S) i— (M, B(M)) is 
called L-measurable if there exists a sequence {/„} of simple function as in Definition [TT| such that 

* fn ^ f a.e. (i.e. c-a.e., where c is the control measure of L), 

• for every A £ S, the sequence of simple integrals {J^fndL} converges in probability, as 

n — >• oo. 

For integrable measurable functions, define 

fdL := P - lim / f^dL. 

A J A 



Rajput & Rosinski ( 1989 1 have pointed out that the above integral is well-defined in the sense 
that it does not depend on the approximating sequence {/n}. Also, the necessary and sufficient 
conditions for the existence of the integral J fdL can be expressed in terms of the characteristics 



of L and can be found in Rajput & Rosinski] ( |1989[ Theorem 2.7), which says the following. Let 



/ : {S, S) !-)• (M, ;B(]R)) be a measurable function. Let L be a Levy basis with CQ (a, b, u{dx, •), c) 
Then / is integrable w.r.t. L if and only if the following three conditions are satisfied: 

/ |Vi(/(z),z)|c(o!z) <oo, / |/(z)|26(z)c((iz) <oo, f V2{f{z),z)c{dz) <oo, (9) 
^ S J S J s 

where for q{x) := xl^_i^^{x), 

F,(u.z):=,.<.(z)+ K.(„,z):= / mi„(l, M>W., z). 

Note that such integrals have been defined for deterministic integrands. However, in the context 
of ambit fields, which we will focus on in this paper, we typically encounter stochastic integrands 
representing stochastic volatility, which tends to be present in most applications we have in mind. 
Since we often work under the independence assumption that the stochastic volatility a and the Levy 
basis L are independent, it has been suggested to work with conditioning to extend the definition by 



Rajput & Rosinski ( 1989 1 to allow for stochastic integrands. An alternative concept, which directly 



allows for stochastic integrands which can be dependent of the Levy basis, is the integration concept 



by Walsh ( 1986 1, which we study next. 



2.2.2 Integration w.r.t. martingale measures introduced by |Walsh| ( |1986^ 



The integration theory due to Walsh (19861 can be regarded as Ito integration extended to random 
fields. In the following we will present the integration theory on a bounded domain and comment 
later on how one can extend the theory to the case of an unbounded domain. 

Here we treat time and space separately, which allows us to work with a natural ordering (intro- 
duced by time) and to relate the integrals w.r.t. to Levy bases to martingale measures. In the following, 
we denote by 5 a bounded Borel set in A' = M'' for a d G No (where d + 1 = k) and S = B{S) 
denotes the Borel cr-algebra on S. Since S is bounded, we have in fact S = B{S) = Bh{S). 
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Let L denote a Levy basis on (5 x [0, T],B{S x [0, T])) for some T > 0. For any A £ Bb(S) 
and < t < T, we define 

Lt[A) = L{A,t) = L{A X (0,t]). 
Here Lt{-) is a measure-valued process, and for a fixed set A G Bi,{S), Lt{A) is an additive process 



in law. In the following, we want to use the Lt{A) as integrators as in Walsh ( 1986 1. In order to do 
that, we work under the square-integrability assumption, i.e.: 

Assumption (Al): For each A £ Bb{S), we have that Lt{A) £ L'^{Q, T, P). 

In the following, we will, unless otherwise stated, work without loss of generality under the zero-mean 
assumption on L, i.e. 

Assumption (A2): For each A £ Bh{S), we have that ¥.{Lt{A)) = 0. 



Next, we define the filtration Ft by 



•0 

t+l/n' 



where T'^ = (j{Ls{A) : A £ Bb{S),0 <s<t}vM, 



(10) 



and where M denotes the P-nuU sets of T. Note that Tt is right-continuous by construction. One 
can show that under the assumptions (Al) and (A2) and for fixed A £ Bb{S), {Lt{A))o<t<T is a 
(square-integrable) martingale with respect to the filtration {Tt)o<t<T- Note that these two properties 
together with the fact that Lq{A) = a.s. ensure that (-^^i(^))t>o AeB(i?<*) i^ ^ martingale measure 
with respect to {Tt)o<t<T in the sense of Walsh ( 1986[ l. Furthermore, we have the following orthogo- 
nality property: If A~B £ Bb{S) with Ar]B = 0, then Lt{A) and Lt{B) are independent. Martingale 
measures which satisfy such an orthogonality property are referred to as orthogonal martingale mea- 
sures by Walsh ( |1986 ), see also Bamdorff-Nielsen, Benth & Veraart| (201 1 1 for more details. Note that 
orthogonal martingale measure are worthy, see Walsh ( I986[ Corollary 2.9), a property which makes 
them suitable as integrators. For such orthogonal martingale measures, [Walsh ( 1986 1 introduces their 
covariance measure Q by 



Q{Ax [0,t]) = <LiA) >t, 



(11) 



for ^ £ B{S). Note that Q is a positive measure and is used by Walsh ( 1986 1 when defining stochastic 
integration with respect to L. 



Walsh (19861 defines stochastic integration in the following way. Let ({S^, s) be an elementary 



random field s), i.e. it has the form 



(12) 



where < o < t, o < 6, X is bounded and Ja-measurable, and A £ S. For such elementary 
functions, the stochastic integral with respect to L can be defined as 



at s) Lidt ds) := X {Ltr^M ^ B) - Lt^A D B)) 



(13) 



JB 



for every B £ S.lt turns out that the stochastic integral becomes a martingale measure itself in B (for 
fixed a, b, A). Clearly, the above integral can easily be generalised to allow for integrands given by 
simple random fields, i.e. finite linear combinations of elementary random fields. Let T denote the set 
of simple random fields and let the predictable cr-algebra V be the cj-algebra generated by T. Then 
we call a random field predictable provided it is "P-measurable. The aim is now to define stochastic 
integrals with respect to L where the integrand is given by a predictable random field. 
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In order to do that Walsh ( 1986 1 defines a norm || • ||l on the predictable random fields ( by 



li 



E 



[0,T]x5 



C\ts)Qidtds) 



(14) 



which determines the Hilbert space Vl '■= L^{Q x [0, T] x 5, V, Q), which is the space of predictable 
random fields C, with ||C|li < oo, and he shows that T is dense in Vl- Hence, in order to define 
the stochastic integral of C G Vl, one can choose an approximating sequence {Cn}n C T such that 
lie — CnlU — as n — )- oo. Clearly, for each S, ^j^^ Cn(^, s) L{d^, ds) is a Cauchy sequence 
in L^(r2, T, P), and thus there exists a limit which is defined as the stochastic integral of C- 

Then, this stochastic integral is again a martingale measure and satisfies the following Ito-type 
isometry: 



E 



^J[0,T]xS 

see ( |Walsh|1986t Theorem 2.5) for more details. 



C{ts)L{dtds) 



li, 



(15) 



2.2.3 Relation between the two integration concepts 



The relation between the two different integration concept has been discussed in Bamdorff-Nielsen, 



Benth & Veraart (201 1 pp. 60-61), hence we only mention it briefly here. 



Note that the Walsh (1986 1 theory defines the stochastic integral as the L^-Umit of simple random 
fields, whereas |Rajput & Rosinski ( [1989| ) work with the P- limit. Bamdorff-Nielsen, Benth & Veraart 



( |2011| l point out that deterministic integrands, which are integrable in the sense of Walsh, are thus also 
integrable in the Rajput and Rosinski sense since the control measure of [Rajput & Rosin^E (1989 1 
and the covariance measure of|Walsh|(|1986[l are equivalent. 



3 General aspects of the theory of ambit fields and processes 

In the following we will show how stochastic processes and random fields can be constructed based on 
Levy bases, which leads us to the general framework of ambit fields. This section reviews the concept 



of ambit fields and ambit processes. For a detailed account on this topic see Bamdorff-Nielsen, Benth 



& Veraart (2011 1 and Barndorff-Nielsen & Schmiegel (2007 1. 



3.1 The general framework 

The general framework for defining an ambit process is as follows. Let Y = {Yt (x)} with lt(x) := 
y (x, t) denoting a stochastic field in space-time A' x M. In most applications, the space X is chosen 
to be M"' for d = 1, 2 or 3. Let w (6) = (x (9) ,t{e)) denote a curve inX xR. The values of the field 
along the curve are then given by = Y^g-^ (x {6)). Clearly, X = {Xg} denotes a stochastic process. 
Further, the stochastic field is assumed to be generated by innovations in space-time with values Yj (x) 
which are supposed to depend only on innovations that occur prior to or at time t and in general only 
on a restricted set of the corresponding part of space-time. I.e., at each point (x, t), the value of Yt (x) 
is only determined by innovations in some subset At (x) of X x Rt (where = (— oo, t]), which 
we call the ambit set associated to (x, t). Furthermore, we refer to Y and X as an ambit field and an 
ambit process, respectively. 

In order to use such general ambit fields in applications, we have to impose some structural as- 
sumptions. More precisely, we will define Yt (x) as a stochastic integral plus a drift term, where the 
integrand in the stochastic integral will consist of a deterministic kernel times a positive random vari- 
ate which is taken to express the volatility of the field Y. More precisely, we think of ambit fields as 
being defined as follows. 
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Definition 13. Using tiie notation introduced above, an ambit field is defined as a random field of the 
form 

Yt{x)=fi + [ h (x, t; s) as (^) L (d^, ds) + /" q (x, t; ^, s) a, (^) d^ds, (16) 

JAi(x) JDt(x) 

provided the integrals exist, where At (x), and Dt (x) are ambit sets, h and q are deterministic func- 
tions, cr > is a stochastic field referred to as volatility or intermittency , a is also a stochastic field, 
and L is a Levy basis. 

Remark 2. Note that compared to the base model ([T]l we introduced in the Introduction, the ambit 
field defined in ([T6]l also comes with a drift term and stochastic volatility introduced in form of a 
stochastic integrand. In Section |5} we will describe in detail how such a stochastic volatility field a 
can be specified and what kind of complementary routes can be taken in order to allow for stochastic 
volatility clustering. 

The corresponding ambit process X along the curve w is then given by 

Xe = fi+ [ h{t{e);ts)as{^)L{dtds)+ [ qit(e);t s)as{Od$ds, (17) 
JA{e) JD(e) 

where A{e) = ^t(,)(x(e)) and D{e) = A(e)(x(e)). 

In Section |3.2[ we will formulate the suitable integrability conditions which guarantee the exis- 
tence of the integrals above. 

Of particular interest in many applications are ambit processes that are stationary in time and 
nonanticipative and homogeneous in space. More specifically, they may be derived from ambit fields 
Y of the form 

Ytix)=f,+ [ h{x-tt-s)as{^)L{dtds)+ [ g (x - ^, t - s) a, (0 d^ds. (18) 

Here the ambit sets At (x) and Dt (x) are taken to be homogeneous and nonanticipative, i.e. At (x) 
is of the form At (x) = ^ + (x, t) where A only involves negative time coordinates, and similarly for 
Dt (x). In addition, a and a are chosen to be stationary in time and space and L to be homogeneous. 



3.2 Integration for general ambit fields 



Ambit fields have initially been defined for deterministic integrands using the Rajput & Rosinski 



( 1989 1 integration concept. Their definition could then be extended to allow for stochastic integrands 
which are independent of the Levy basis by a conditioning argument. As discussed before, the integra- 
tion framework developed by Walsh (1986]) has the advantage that it allows for stochastic integrands 
which are potentially dependent of the Levy basis and enables us to study dynamic properties (such 
as martingale properties). Let us explain in more detail how the Walsh ( |1986| ) integration concept can 
be used to define ambit fields using an Ito-type integration concept. 



One concern regarding the applicability of the Walsh ( 1986 1 framework to ambit fields might be 
that general ambit sets At{x) are not necessarily bounded, and we have only presented the Walsh] 
(19861 concept for a bounded domain. However, the stochastic integration concept reviewed above 



can be extended to unbounded ambit sets using standard arguments, cf. Walsh ( 1986 p. 289). Also, 
as pointed out in Walsh ( I986[ p. 292 ), it is possible to extend the Walsh ( 19861 integration concept 
beyond the L^-framework, cf. Walsh] (] 1986] p. 292). 



Note that the classical Walsh (1986 1 framework works under the zero mean assumption, which 
might not be satisfied for general ambit fields. However, we can always define a new Levy ba sis L by 
setting L := L — E(L), which clearly has zero mean. Then we can define the Walsh ( 1986) integral 



11 



3 GENERAL ASPECTS OF THE THEORY OF AMBIT FIELDS AND PROCESSES 



w.r.t. L, and we obtain an additional drift term which needs to satisfy an additional integrability 
condition. 

However, the main point we need to address is the fact that the integrand in the ambit field does 
not seem to comply with the structure of the integrand in the Walsh-theory. More precisely, for ambit 
fields with ambit sets At(x) C M.'^ x (— oo,t], we would like to define Walsh-type integrals for 
integrands of the form 

ats) :=C(x,t;^,s) :=lA,i^){ts)h{^,t;ts)aM. (19) 
The original Walsh's integration theory covers integrands which do not depend on the time index t. 



Clearly, the integrand given in ( 19 1 generally exhibits t-dependence due to the choice of the ambit set 
At (x) and due to the deterministic kernel function h. 

Suppose we are in the simple case where the ambit set can be represented as At{x.) = Bx (— oo, t], 
where B G B{M.'^) does not depend on t, and where the kernel function does not depend on t, i.e. 
/i(x, t; ^, s) = /i(x; ^, s). Then the Walsh-theory is directly applicable, and provided the integrand is 
indeed Walsh-integrable, then (for fixed B and fixed x) the process 



h{x; ^, s)asi^)L{d^, ds) 



B 

is a martingale. 

Note that the i-dependence (and also the additional x-dependence) for general integrands in the 
ambit field is in the deterministic part of the integrand only, i.e. in IlAt(x)(^5 s)/i(x, t; ^, s). Now in 
order to allow for time t- (and x-) dependence in the integrand, we can define the integrals in the 
Walsh sense for any fixed t and for fixed x. Note that we treat x as an additional parameter which 
does not have an influence on the structural properties of the integral as a stochastic process in t. 

It is clear that in the case of having t-dependence in the integrand, the resulting stochastic integral 
is, in general, not a martingale measure any more. However, the properties of adaptedness, square- 
integrability and countable additivity carry over to the process 



({x,t;ts)L{dtds) 

(for fixed x) since it is the L^-limit of a stochastic process with the above mentioned properties. 



In order to ensure that the ambit fields (as defined in ( 16 1) are well-defined (in the Walsh-sense), 



throughout the rest of the paper we will work under the following assumption: 

Assumption 1. Let L denote a Levy basis on Sx (— oo, T], where S denotes a not necessarily bounded 
BorelsetSinX = for some d £ IN. Define the new Levy basis L := L — E(L). We extend the 



definition of the covariance measure Q of L, see {11 \ to an unbounded domain and, next, we define 



a Hilbert space Vj^ with norm \ \ ■ \ as in {14 \ (extended to an unbounded domain) and, hence, we 
have an ltd isometry of type (15) extended to an unbounded domain. We assume that, for fixed x and 
t, 

satisfies 

1. C G Vj^, 



2. \\C\\j^2=E 



^- k>,x ds) <oo 



< oo. 
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Remark 3. Note that alternatively, we could work with the cadlag elementary random fields 



where X{u}) is assumed to be J^^-adapted and the remaining notation is as in ( 12 1. Next, one can 
construct a cr-algebra from the corresponding simple random fields and one would then define the 
stochastic integral for C,*{^, s—,uj), since clearly adaptedness and the cadlag property of C,*{^,s,uj) 
implies predictability of C* s— , cj). 

3.3 Cumulant function for stochastic integrals w.r.t. a Levy basis 

Next we study some of the fundamental properties of ambit fields. Throughout this subsection, we 
work under the following assumption: 

Assumption 2. The stochastic fields a and a are independent of the Levy basis L. 

Now we have all the tools at hand which are needed to compute the conditional characteristic func- 
tion of ambit fields defined in ( [T6| ) where a and L are assumed independent and where we condition 
on the cr-algebra F'^ = -T^f (x) which is generated by the history of a, i.e. 

JT(x)=aK(0:(^,5)GAt(x)}. 

Proposition 1. Assume that Assumption^holds. Let C" denote the conditional cumulant function 
when we condition on the volatility field a. The conditional cumulant function of the ambit field 



defined by {16) is given by 



CUet [ h{^,t;ts)as{^)L{dCds)\ 

= log l^exp {^^ f^^^ ^(^' ^' s)as{$)L{di, ds)^ F''^ ^ (20) 
C {0/i(x, t; t s)asiO t L'it s)} cidt ds), 



where L' denotes the Levy seed and c is the control measure associated with the Levy basis L, cf. 
and 



Proof. The proof of the Proposition is an immediate consequence of [Rajput & RosinsE (1989] Propo 



sition 2.6). □ 



Corollary 1. In the case where L is a homogeneous Levy basis, equation {20) simplifies to 



C'lexf h{^,t;ts)asU)Lidtds)\ = [ C {eh{x,t; ^ s)as{^) t L'} d^ds. 
[ JAt{^) J JAt(x) 

3.4 Second order structure of ambit fields 

Next we study the second order structure of ambit fields. Throughout the Section, let 

Yt{ic)= [ h{^,t;^,s)a,{OL{dtds), (21) 

where a is independent of L, i.e. Assumption [IJholds, and L' is the Levy seed associated with L. 
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Proposition 2. Let t, t, x, x > and let lt(x) be an ambit field as defined in {21 ). The second order 
structure is then as follows. The means are given by 



¥.{Yt{^)\Fn= I h{^,t;^,s)as{m{L'{ts))c{dtds), 

JAt(x) 

E (yt(x)) = / /i(x, t; t s)E {asm HL'{t sMd^ ds). 
The variances are given by 

Var{Yt{^)\F'')= [ h\^,t; ^ s)a^MVar{L' {t s))c{dt ds), 

Var (yt(x)) = / /i2(x, t- 1 s)E {a^,{^)) Var{L'{^, s))c{dt ds) 
JAt{yi) 

+ I [ hi^,t;ts)hi^,t;lrs)p{srs,tl)nL\ts))EiL\ls))c{dtds)cidld~s), 

JAtix.) JAtix.) 

where p{s, s, ^, ^) = E (^crs{$)o's{(,)^ — E (o"s(^)) E The covariances are given by 

Cov {Yt{^),Yi{i^)\Fn = [ h{^,t;ts)h{±,i;^,s)a^^{^)Var{L'{ts))c{dtds), 

JAt{x)r\Ai{k) 

Cov{Yt{^),Yi{i^)) = [ h{^,t;ts)h{5,,i;ts)E{al{^))Var{L'{ts))c{dtds) 

JAt(x)nAf(x) 

+ f f h{^,t;ts)h{5^,i;ls)p{s,s,^,m{L'{ts))E{L'{l~s))c{dld~s)c{dCds). 

JAt(x) JAr(x) 



Proof. From the conditional cumulant function (20 1, we can easily deduce the second order structure 



conditional on the stochastic volatility. Integrating over a and using the law of total variance and 
covariance leads to the corresponding unconditional results. □ 

Note that it is straightforward to generalise the above results to allow for an additional drift term 



as in ( 16 1. 



The second order structure provides us with some valuable insight into the autocorrelation struc- 
ture of an ambit field. Knowledge of the autocorrelation structure can help us to study smoothness 
properties of an ambit field, as we do in the following section. Also, from a more practical point of 
view, we could think of specifying a fully parametric model based on the ambit field. Then the sec- 
ond order structure could be used e.g. in a (quasi-) maximum-likelihood set-up to estimate the model 
parameters. 

3.5 Smoothness conditions 

Let us study sufficient conditions which ensure smoothness of an ambit field. 
3.5.1 Some related results in the literature 



In the purely temporal (or null-spatial) case, which we will discuss in more detail in Section 3.7 



smoothness conditions for so-called Volterra processes have been studied before. In particular, De 



creusefond| ( j2002) shows that under mild integrability assumptions on a progressively measurable 



stochastic volatility process, the sample-paths of the volatility modulated Brownian-driven Volterra 



process are a.s. Holder-continuous even for some singular deterministic kernels. Note that Decreuse 



fond ( 2002 1 does not use the term stochastic volatility in his article, but the stochastic integrand he 
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considers could be regarded as a stochastic volatility process. Also, Mytnik & Neuman (2011 1 study 
sample path properties of Volterra processes. 

In the tempo-spatial context, or generally for random fields, smoothness conditions have been 



discussed in detail in the literature. For textbook treatments see e.g. Adler ( 1981 1, Adler & Taylor 



(20071 and Azais & Wschebor ( 2009[ ). Important articles in this context include the following ones. 
Kent (19891 formulates sufficient conditions on the covariance function of a stationary real- valued 
random field which ensure sample path continuity. 



Rosinski (1989 1 studies the relationship between the sample-path properties of an infinitely di- 
visible integral process and the properties of the sections of the deterministic kernel. The study is 
carried out under the assumption of absence of a Gaussian component. In particular, he shows that 
various properties of the section are inherited by the paths of the process, which include boundedness, 
continuity, differentiabihty, integrability and boundedness of pth variation. 



Marcus & Rosinskil ( 2005 ) extend the previous results to derive sufficient conditions for bound- 



edness and continuity for stochastically continuous infinitely divisible processes, without Gaussian 
component. 



3.5.2 Sufficient condition on the covariance function 

In the following, we write 

p{t, x; i, x) = Cov (yt(x), y^(i)) , 



where the covariance is given as in Proposition [2] We can apply the key results derived in Kent ( 1989 1 
to ambit fields. 

Let t,hi € M and x, h2 € and h := (/ii,h2)'. For each (x, f)' we assume that p{t + 
/ii, X + h2; t — /ii, X — h2) is A;-times continuously differentiable with respect to h for a A; G INq. 
We write pfc(h; t, x) for the polynomial of degree k which is obtained from a Taylor expansion of 

+ /ii, X + h2; t — /ii, X — 112) about h = for each (x, t). In the following, we denote by || • || 
the Euclidean norm. 

Propositions. Foreach {x,ty suppose that p{t+hi,x+h2;t— hi, x— 112) is d+l-times continuously 
differentiable with respect to h and that there exists a constant 7 > such that 



( l|h||'^+^ 

sup{|p(t + /ii,x + h2;t-/ii,x-h2) - pd+i(h; t, x)|} = O m, .mo ,^ 

(x,t) Vliog||h|ir+^ 



(22) 



1 1 h| I — )• 0, where the supremum is computed over (x, t) being in each compact subset o/R'^+^. Then 
there exists a version of the random field {l((x), (x, t) G M'^+^j which has almost surely continuous 
sample paths. 



Proof. The result is a direct consequence of Kent ( 1989 Theorem 1 and Remark 6). 



□ 



Remark 4. As pointed out in Kent ( 1989 Remark 3), ( |22| ) could be replaced by the stronger con- 
ditions of the supremum being of order 0(||h||'^+^+'') as h — for some constant /3 > 0. This 



condition is easier to check in practice and if it holds for any /3 > 0, then it implies that ( [22] ) holds for 
all 7 > 0. 



Remark 5. As pointed out in Kent ( 1989 Remark 5) and Adler ( 1981 p. 60), as soon as we have a 
Gaussian field, milder conditions ensure continuity. Note that an ambit field is Gaussian if the Levy 
basis is Gaussian and the stochastic volatility component is absent (or purely deterministic). 
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3.6 Semimartingale conditions 

Next we derive sufficient conditions which ensure that an ambit field is a semimartingale in time. 
This is interesting since in financial applications we typically want to stay within the semimartingale 
framework whereas in applications to turbulence one typically focuses on non-semimartingales. 

We will see that a sufficient condition for a semimartingale is linked to a smoothness condition 
on the kernel function. When studying the semimartingale condition, we focus on ambit sets which 
factorise as ^t(x) = [0,t] x A{x.), which is in hne with the [Walsh, 1 986| )-framework. We start with 
a preliminary Lemma. 

Lemma 1. Let L be a Levy basis satisfying (AI ) and (A2) and a be a predictable stochastic volatility 
field which is integrable w.rt. L. Then 

Mt{A{-^))= f f as{$)L{dtds), (23) 

Jo JAix) 

is an orthogonal martingale measure. 

Proof. See Walsh (1986} p. 296) for a proof of the lemma above. □ 



Assumption 3. Assume that the deterministic function u i— )• n, •, •) is differentiable (in the second 
component) and denote by h' = •, •) the derivative with respect to the second component. 

Then for s <twe have the representation 



h{x,t;i,,s) = h{x,s;$,s) + J h'{x,u;i,, s)du, 



(24) 



pro vided that h(x , s; ^, s) exists. Further, we assume that both components in the representation ( (24 



are 



Walsh 



(1986]-integrable w.rt. M and Ci(^5 s) = /i(x, s; ^, s) and C2(^7 s) = h'{x, u; ^, s)d 



satisfy Ci e Vm, WdWu^ < oo, L^.^ \Ci{i,s)\&M{d$„ds) < oo fori = 1,2. 



Proposition 4. Let M be defined as in ( |23| ) with covariance measure Q m, ond assume that h satisfies 
Assumption ^and 



E / \F {u, x] s, $,{)F {u,yi; s, i2)\Q M{dii, d$,2, ds)du < oo, 

J SxSx[0,T\xG 

for 

F{u,x;s,$,) = h'{x,u;^,s)I[s,t]{u)- 
Then (J^ /^^^^ /i(x, t; ^, s)AI{d^, ds))t>o ^ semimartingale with representation 



h{t,x;s,$)M{d^,ds) 
JA{-s.) 

= 11 h{s,x;s,^)M{d$„ds) + / / / x; s, ^)M((i^, ds)dii. (25) 

Jo J A{-x.) Jo Jo Ja{x.) 



Clearly, the first term in representation ( [251 ) is a martingale measure in the sense of Walsh and the 
second term is a finite variation process. 



Proof. The result follows from the stochastic Fubini theorem for martingale measures, see Walsh 



( [1986 Theorem 2.6). In the following, we check that the conditions of the stochastic Fubini theo- 



rem are satisfied. Let {G,G, Leb) denote a finite measure space. Concretely, take G = [0,T] and 
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Q = B{G). Note that u i— )• h{-,u, •, •) has finite first derivative at all points, hence its derivative is 
^-measurable. Also, the indicator function is ^-measurable since the corresponding interval is an 
element of Q. Overall we have that the function 



is ^ X P-measurable and 



F(n,x;s,^) = h' {^,u; s)\,t]{u) 



E / \F{u, x;s, $i)F{u, :>c;s, ^2)\QM{d$i, d$2, ds)du < oo, 

JsxSx[0,T]xG 

where Qm is the covariance measure of M. Then 



h{t,x;s,^)M{d$,ds) 
JA{^) 

[ h{s,:xi;s,$)M{d$„ds)+ [ [ [ h'{u,:x.;s,$)duM{d$„ds) 

Ja(x.) Jo Ja{x.) Js 

t f ft ru p 

I h{s,x;s,$,)M{d$„ds) + / / h'{u,x; s,$)M{d$„ds)du. 

JA(x) Jo Jo J A{yi) 



□ 



3.7 The purely-temporal case: Volatility modulated Volterra processes 

The purely temporal, i.e. the null-spatial, case of an ambit field has been studied in detail in recent 
years. Here we denote by L = {Lt)teR a Levy process on M. Then, the null-spatial ambit field is in 
fact a volatility modulated Levy-driven Volterra (VAiCV) process denoted by y = {Yt}^^^, where 



q{t,s)asds, 



(26) 



where /i is a constant, k and q are are real- valued measurable function on M^, such that the integrals 
above exist with k {t, s) = q {t, s) = for s > t, and a and a are cadlag processes. 

Of particular interest, are typically semi-stationary processes, i.e. the case when the kernel func- 
tion depends on t and s only through the difference t — s. This determines the class of Levy semista- 
tionary processes (CSS), see Bamdorff-Nielsen, Benth & Veraart ( 2012[ ). Specifically, 



Yt = fj.+ k(t- s)as-dZs + / q{t - s)asds, 



(27) 



k and q are non-negative deterministic functions on M with k (x) = q{x) = for x < 0. Note 
that an CSS process is stationary as soon as a and a are stationary processes. In the case that L is 



a Brownian motion, we call Y a Brownian semistationary (BSS) process, see Bamdorff-Nielsen & 



Schmiegel| ( |2009l ). 

The class of BSS processes has been used by Bamdorff-Nielsen & Schmiegel| ( [2009| l to model tur- 
bulence in physics. In that context, intermittency, which is modelled by a, plays a key role, which has 
triggered detailed research on the question of how intermittency can be estimated non-parametrically. 



Recent research, see e.g. |Bamdorff-Nielsen et al. (20091, Bamdorff-Nielsen, Corcuera & Podolskij 



(2011 2012 1, has developed realised multipower variation and related concepts to tackle this impor- 
tant question. 

The class of CSS processes has subsequently been found to be suitable for modelling energy 
spot prices, see Bamdorff-Nielsen, Benth & Veraart ( 2012| |, [Veraart & Veraart (2012 1. Moreover, 



Bamdorff-Nielsen, Benth, Pedersen & Veraart| ( |2012[ ) have recently developed an anticipative stochas- 
tic integration theory with respect to VM.CV processes. 
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4 Illustrative example: Trawl processes 

Ambit fields and processes constitute a very flexible class for modelling a variety of tempo-spatial 
phenomena. This Section will focus on one particular sub-class of ambit processes, which has recently 
been used in application in turbulence and finance. 

4.1 Definition and general properties 

Trawl processes are stochastic processes defined in terms of tempo-spatial Levy bases. They have re- 



cently been introduced in Bamdorff-Nielsen ( 201 1 1 as a class of stationary infinitely divisible stochas- 
tic processes. 

Definition 14. Let L be a homogeneous Levy basis on M"^ X M for d G IN. Then, using the same 
notation as before, 

Bb{R'^ X M) := {A e B{R'^ X M) : Leb{A) < oo}. 
Further, for an yl = G x M), we define At:= A + (0, t). Then 

Yt= [ lAits-t)L{dtds) = L{At), (28) 

defines the trawl process associated with the Levy basis L and the trawl A. 



The assumption Leb{A) < oo in the definition above ensures the existence of the integral in (28 1 
as defined by [Rajput & Rosin^E ( 1989 1. 



Remark 6. If A C (— oo, 0] x M'', then the trawl process belongs to the class of ambit processes. 

The intuition and also the name of the trawl process comes from the idea that we have a certain 
tempo-spatial set - the trawl - which is relevant for our object of interest. I.e. the object of interest at 
time t is modelled as the Levy basis evaluated over the trawl At. As time t progresses, we pull along 
the trawl (like a fishing net) and hence obtain a stochastic process (in time t). For the time being, we 
have in mind that the shape of the trawl does not change as time progresses, i.e. that the process is 
stationary. This assumption can be relaxed as we will discuss in Section [5] 

Example 7. Let d = 1 and suppose that the trawl is given by At = {{x,s) : s < t,0 < x < 
exp(— 0.7(t — s))}. Figure^illustrates the basic framework for such a process. It depicts the trawl 
at different times t G {2, 5}. The value of the process at time t is then determined by evaluating the 
corresponding Levy basis over the trawl At. 

4.2 Cumulants and correlation structure 



From Bamdorff-Nielsen ( 2011| ) we know that the trawl process defined above is a strictly stationary 



stochastic process, and we can easily derive the cumulant transform of a trawl process, which is given 
by 

C7(C X Yt) := E(exp(iC>t)) = Leb{A)C{C t L') := Xe6(A)E(exp(iCi'))- (29) 



From the equation ( [29] ), we see immediately that the law of the trawl process is infinitely divisible 
since the corresponding Levy seed has infinitely divisible law. 

Remark 7. To any infinitely divisible law vr there exists a stationary trawl process having tt as its 



one-dimensional marginal law. This follows from formula (29 1, cf. Bamdorff-Nielsen (2011 1 
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Figure 1: Example of a relevant choice of the trawl: At = {{x, s) : s < t,0 < x < exp(— 0.7(t— s))}. 
The shape of the trawl does not change as time t progresses and, hence, we obtain a stationary process. 
The value of the process is obtained by evaluating L{At) for each t. 



We can now easily derive the cumulants of the trawl process which are given by Ki{Yt) = 
Leb{A)ni{L'), for i G IN, provided they exist. In particular, the mean and variance are given by 

E(yt) = Leb{A) E(L'), Var{Yt) = Leb{A) Var{L'). 

Marginally, we see that the precise shape of the trawl does not have any impact on the distribution 
of the process. The quantity which matters here is the size, i.e. the Lebesgue measure, of the trawl. So 
two different specifications of the ambit set, which have the same size, are not identified based on the 
marginal distribution only. 

However, we will see that the shape of the trawl determines the autocorrelation function. More 
precisely, the autocorrelation structure is given as follows. Let h > 0, then 

p{h) := Cov{Yt,Yt+h) = Leh{A n Ah)Var{L'). (30) 

For the autocorrelation, we get 

ih\ n (V V \ Leb{Ar\Ah) 
r{h)=Cor{Yt,Yt^,) = ^-^. 



4.3 Levy-Ito decomposition for trawl processes 

In the following we study some of the sample path properties of trawl processes. First of all, we study 
a representation result for trawl process Y, where we split the process into a drift part, a Gaussian 
part and a jump part in a similar fashion as in the classical Levy-Ito decomposition. Recall that L 
is a homogeneous Levy basis on X M. In the following, s denotes the one dimensional temporal 
variable and ^ denotes the d-dimensional spatial variable. 

From the Levy-Ito decomposition, see Q, we get the following representation result for the trawl 
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process Yt = L{At) defined in p8| ), i.e 



L{At) = aLeb{A) + W{At)+ [ xN{dx,At)+[ x{N - n){dx, At) 

J{\x\>l} J{\x\<i} 

= aLeb{A) + W{At)+ [ [ xN{dx,d$,ds) + [ [ x{N - n){dx,d$,ds), 

J At ^{|x|>l} J At J{\x\<l} 

(31) 

where W{At) ~ iV(0, bLeb{A)) and n{dx, d^, ds) = v[dx)d$,ds and Leb{At) = Leb{A). 

Example 8. Suppose the trawl process Yt is defined based on a Levy basis with characteristic quadru- 
plet (0, b, 0, Leb). Then it can be written as 

Yt = W{At)r^ N{Q,bLeb{A)). 

Assume fiirther that d = 1 and that the trawl is given by At = {{x, s) : s < t,0 < x < exp(— A(i — 
s))}, for a positive constant A > 0. Then Leb{At) = j, hence Yt = W{At) ~ N(0, j) and the 
autocorrelation function is given by r(h) = exp(— A/i), for h > 0. 

4.4 Generalised cumulant functional 

We have already studied the cumulant function of Levy bases and ambit fields. Here, we will in 
addition focus on the more general cumulant functional of a trawl process, which sheds some light on 
important properties of trawl processes. 

Definition 15. Let Y = (Yt) denote a stochastic process and let /j, denote any non-random measure 
such that 



fi{Y) = / Ytf^idt) < oo, 

JR 

where the integral should exist a.s.. The generalised cumulant functional of Y w.r.t. fi is defined as 

C{0tMy)} = logE(exp {iefi{Y))). 
When we compute the cumulant functional for a trawl process, we obtain the following result. 

Proposition 5. Let Yt = L{At) denote a trawl process and let fi denote any non-random measure 
such that /u(y) = J^Yt fi{dt) < oo, a. s.. Given the trawl A, we will further assume that for all 



hAit s) = (^j^ 1a(^, s - t)fiidt)^ < oo, 



and that hA{$,, s) is integrable with respect to the Levy basis L. Then the cumulant function of fj,(Y) 
is given by 

C{etf^iY)} = iea [ hA{i,s)dids-\e% f h\{$,s)d$ds 

+ / / [exp{i9ux) — 1 — 1] (x)) i'{dx)x{du), 



XII 



where x i^ the measure on M obtained by lifting the Lebesgue measure onW^ xM^toM. by the mapping 

{i,s)^{i,hA{s,i)). 
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Proof. An application of Fubini's theorem (see e.g. Bamdorff -Nielsen & Basse-O'Connor (2011 1) 
yields 



[ [ lA{^,s-t)Lid^,ds)fi{dt)= [ (I lA{^,s-t)fi{dt)] L{d^^ds) 

^ V ' 

= [ hAits)L{dtds). 

From this we find that the cumulant function of fJ,{Y), i.e. the generalised cumulant functional of Y 
w.r.t. /i, is given by 

c{etKy)}= I c{9hA{i,s)XL']dids 

= iOa [ hA{^, s)d$ds - h^b [ h\{$, s)d$ds 

JrxR'^ 2 jRxR'i 

+ / [ex.p{i9hA{^, s)x) — 1 — i9hA{^, s)xl[^i i]{x)) u{dx)d$ds. 

JMxRd Jr 

Note that the latter part, i.e. the jump part of C{0 | fi{Y)}, can be recast as 

/ / (exp{i9hA{^, s)x) — 1 — i9hA{$, v{dx)d$,ds 

jRxR'i Jr 

[ex.-p{i9ux) — 1 — z0uxl[_i v{dx)x{du), 



XJl 



where x is the measure defined as above. Then the result follows. □ 

In the following, we mention three relevant choices of the measure ^. Let 5t denote the Dirac 
measure at t. We start with a very simple case: 

Example 9. Suppose fi{t) = 6 (dt) for a fixed t. Hence 

hA{i, s) = (^j^ lA{t s - t)fi{dt)^ = 1a{^, s - t) < oo, 

and clearly, Ha^^, s) is integrable with respect to L. Then 

C{9 t n{Y)} = C{9 t Yt} = Leb{A)C{9 t L'}. 



This is exactly the result we derived in Section |4.2| above. More interesting is the case when n 
is given by a linear combination of different Dirac measures, since this allows us to derive the joint 
finite dimensional laws of the trawl process and not just the distribution for fixed t. 

Example 10. Suppose /i {dt) = 9i6t^ (dt) + • • • + 9ndt„ (dt) for constants 9i, . . . ,9n G M and 
times ti < ■ ■ ■ < tn for n G JN. As in the example before, the integrability conditions which were 



needed to derive (32 ) are satisfied and, hence, ( |j2| ) gives us the cumulant function of the joint law of 
Finally, another case of interest is the integrated trawl process, which we study in the next example. 



Example 11. Let fj, (dt) = 1/ (t) dtfor an interval I o/M. Then {32 \ determines the law of jj Ygds. 



Remark 8. The last example is particularly relevant if the trawl process is for instance used for 
modelling stochastic volatility. Note that such an application is feasible since the trawl process is 
stationary and we can formulate assumptions which would ensure the positivity of the process as well 
(e.g. if we work with a Levy subordinator as the corresponding Levy seed). In that context, integrated 
volatility is a quantity of key interest. 
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4.5 The increment process 

Finally, we focus on the increments of a trawl process. Note that whatever the type of trawl, we have 
the following representation for the increments of the process for s < t, 

Yt-Ys = L {At\As) - L {As\At) , almost surely. (33) 

Due to the independence of L {At\As) and —L {As\At), we get the following representation for the 
cumulant function of the returns 

C{C t Yt - Ys) = C{CtL {At\As)) + C{-C t L {As\At)) 

= Leb{At\As)C{C t L') + LebiAMt)Ci-C t L') (34) 
= Leh[At-s\Ao)C{C t L') + Leb{A^\At.s)C{-(, \ L'). 



4.6 Applications of trawl processes 

Trawl processes constitute a class of stationary infinitely divisible stochastic processes and can be 
used in various applications. E.g. we already pointed out above that they could be used for modelling 



stochastic volatility or intermittency. In a recent article by Bamdorff-Nielsen, Lunde, Shephard & Ve- 



raart 



( 2012| ), integer- valued trawl processes have been used to model count data or integer- valued data 
which are serially dependent. Speaking generally, trawl processes can be viewed as a flexible class 
of stochastic processes which can be used to model stationary time series data, where the marginal 
distribution and the autocorrelation structure can be modelled independently from each other. 



5 Tempo-spatial stochastic volatility/intermittency 

Stochastic volatility/intermittency plays a key role in various applications including turbulence and 
finance. While a variety of purely temporal stochastic volatility/intermittency models can be found in 
the literature, suitable tempo-spatial stochastic volatility/intermittency models still need to be devel- 
oped. 

Volatility modulation within the framework of an ambit field can be achieved by four comple- 
mentary methods: By introducing a stochastic integrand (the term (Tt(x) in the definition of an ambit 
field), or by (extended) subordination or by probability mixing or Levy mixing. We will discuss all 
four methods in the following. 

5.1 Volatility modulation via a stochastic integrand 

Stochastic volatility in form of a stochastic integrand has already been included in the initial definition 
of an ambit field, see ([16]). The interesting aspect, which we have not addressed yet, is how a model 
for the stochastic volatility field at{x) can be specified in practice. We will discuss several relevant 
choices in more detail in the following. 

There are essentially two approaches which can be used for constructing a relevant stochastic 
volatility field: Either one specifies the stochastic volatiUty field directly as a random field (e.g. as 
another ambit field), or one starts from a purely temporal (or spatial) stochastic volatility process and 
then generalises the stochastic process to a random field in a suitable way. In the following, we will 
present examples of both types of construction. 

5.1.1 Kernel-smoothing of a Levy basis 

First, we focus on the modelling approach where we directly specify a random field for the volatility 
field. A natural starting point for modelling the volatility is given by kernel-smoothing of a homo- 
geneous Levy basis - possibly combined with a (nonlinear) transformation to ensure positivity. For 
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instance, let 

a2(x) = v( [ j{^,tt- s)L''{di, ds)\ , (35) 

where L"^ is a homogeneous Levy basis independent of L, j : ^ is an integrable kernel 

function satisfying j(x, u, ^) = for u < and y : M — M_|_ is a continuous, non-negative function. 



Note that cr^ defined by (35 1 is stationary in the temporal dimension. As soon as j(x, — s) 



j*(x — $,t — s) for some function j* in ([35]l, then the stochastic volatility is both stationary in time 
and homogeneous in space. 

Clearly, the kernel function j determines the tempo-spatial autocorrelation structure of the volatil- 
ity field. 

Let us discuss some examples next. 

Example 12 (Tempo-spatial trawl processes). Suppose the kernel function is given by 

j(x,^,t-s) = Ia-(^-x, s-t), 

where C M'^ x (— cxd, 0]. Further, let (x) = A'^ + (x, t); hence ^t'(x) is a homogeneous and 
nonanticipative ambit set. Then 

aU^) = v( [ rA^-^,s-t)L^idtds)) = vl [ L'^idtds)] = V (L^ {A^ {^))) . 

Note that the random field L'^(A^(x)) can be regarded as a tempo-spatial trawl process. 
Example 13. Let 

jXx, ^, t - s) = j*(x, t - s)Ia<t(x)(^, s - t), 

for an integrable kernel function j* and where A'^(x) C M"^ x (— oo,0]. Further, let A1{'x) = 
A''{yi) + (0,t). Then 

a2(x) = F ( / j-(^,x,t-s)L-(dC,(is) I , (36) 

which is a transformation of an ambit field (without stochastic volatility). 

Let us look at some more concrete examples for the stochastic volatility field. 

Example 14. A rather simple specification is given by choosing L'^ to be a standard normal Levy 
basis and V{x) = x^. Then o"^ (^) would be positive and pointwise -distributed with one degree of 
freedom. 

Example 15. One could also work with a general Levy basis, in particular Gaussian, and V given by 
the exponential function, see e.g. Barndorff-Nielsen & Schmiegel\ ( 2004 ) and Schmiegel et al. ( 2QQ5| ). 



Example 16. A non-Gaussian example would be to choose L'^ as an inverse Gaussian Levy basis and 
V to be the identity function. 

Example 17. We have already mentioned that the kernel function j determines the autocorrelation 
structure of the volatility field. E.g. in the absence of spatial correlation one could start off with the 
choice j{x,t — s,^) = exp{— X{t — s)) for X > mimicking the Omstein-Uhlenbeck-based stochastic 
volatility models, see e.g. Barndorff-Nielsen & Schmiegel l\2004^ . 
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5,1.2 Ornstein-Uhlenbeck volatility fields 

Next, we show how to construct a stochastic volatihty field by extending a stochastic process by a 
spatial dimension. Note that our objective is to construct a stochastic volatility field which is stationary 
(at least in the temporal direction). Clearly, there are many possibilities on how this can be done and 
we focus on a particularly relevant one in the following, namely the Omstein-Uhlenbeck-type volatility 
field (OUTVF). The choice of using an OU process as the stationary base component is motivated by 



the fact that non-Gaussian OU-based stochastic volatility models, as e.g. studied in Bamdorff-Nielsen 



& Shephard| ( |2001 1, are analytically tractable and tend to perform well in practice, at least in the 



purely temporal case. In the following, we will restrict our attention to the case d = \, i.e. that the 
spatial dimension is one-dimensional. 

Suppose now that y is a positive OU type process with rate parameter A > and generated by a 
Levy subordinator Y, i.e. 



Yt= f e-^(*-^)dy„ 



We call a stochastic volatility field cr| (x) on M X M an Omstein-Uhlenbeck-type volatility field 
(OUTVF), if it is defined as follows 



n (x) = (x) = e-^'-'Yt + e-'^^^-^^dZ^,, 

Jo 



(37) 



where /i > is the spatial rate parameter and where Z = {Z.^f^^^^ is a family of Levy processes, 
which we define more precisely in the next but one paragraph. 

Note that in the above construction, we start from an OU process in time. In particular, Tt{0) is 
an OU process. The spatial structure is then introduced by two components: First, we we add an 
exponential weight e~^^ in the spatial direction, which reaches its maximum for x = and decays 
the further away we get from the purely temporal case. Second, an integral is added which resembles 
an OU-type process in the spatial variable x. However, note here that the integration starts from 
rather than from — oo, and hence the resulting component is not stationary in the spatial variable x. 
(This could be changed if required in a particular application.) 

Let us now focus in more detail on how to define the family of Levy processes Z. Suppose 

X = < > is a stationary, positive and infinitely divisible process on M. Next we define Z\. = 

as the so-called Levy supra-process generated by X, that is is a family of 

stationary processes such that Z\. has independent increments, i.e. for any < xi < X2 < • • • < x„ 
the processes Z^^^., Z^^l- ~ ^xi\-^ • • • > ^x„\- ~ ^x„_^\. are mutually independent, and such that for each 
X the cumulant functional of equals x times the cumulant functional of X, i.e. 

C{mt Z^\.} = xC{mtX}, 

where 

C{mtX} = logE{e*'"(^)}, 

with m (^X^ = J XsTTi (ds), m denoting an 'arbitrary' signed measure on M. Then at any t G M the 
values Z^^f of Z.^. at time t as x runs through M+ constitute a Levy process that we denote by 



This is the Levy process occurring in the integral in (37 1. 

Note that r is stationary in t and that (x) — )• as x — )■ 0. 

Example 18. Now suppose, for simplicity, that X is an OU process with rate parameter k and gen- 
erated by a Levy process X. Then 



1 



Covin (x),Tt> (x')} = - Var{yi}A-ie-^^l*-* 1^-^^^+^ ^ + Var{Xo}/i 



-Var{Xo}/i~'e 



1 ~K\t-t'\-^i{x+x'y 
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If, furthermore, Var{Yi} = VarjXo} and k = X = fj, then for fixed x and x' the autocorrelation 
function of t is 

Co^{Tt{x),Tt> {x')] = e-''l*-*'le^'^l^-^'l. 

This type of construction can of course be generalised in a variety of ways, including dependence 
between X and Y and also superposition of OU processes. 

Note that the process Tt{x) is in general not predictable, which is disadvantageous given that we 
want to construct Walsh-type stochastic integrals. However, if we choose X to be of OU type, then 
we obtain a predictable stochastic volatihty process. 



5.2 Extended subordination and meta-times 



An alternative way of volatility modulation is by means of (extended) subordination. Extended sub- 
ordination and meta-times are important concepts in the ambit framework, which have recently been 
introduced by Barndorff -Nielsen (2010) and Bamdorff-Nielsen & Pedersen (2012), and we will re- 
view their main results in the following. Note that extended subordination generalises the classical 
concept of subordination of Levy processes to subordination of Levy bases. This in turn will be based 
on a concept of meta-times. 



5.2.1 Meta-times 

This section reviews the concept of meta-times, which we will link to the idea of extended subordina- 
tion in the following section. 

Definition 16. Let S" be a Borel set in M'^. A meta-time T on is a mapping from Bh{S) into Bb{S) 
such that 

1. T {A) and T (B) are disjoint whenever A,B£ Bb{S) are disjoint. 

2. T (u^^i^j) = U°^iT (Aj) whenever Ai, Aa, • • • G BbiS) are disjoint and U^^^^j G BbiS). 
A slightly more general definition is the following one. 

Definition 17. Let be a Borel set in M^. A full meta-time T on 5 is a mapping from B{S) into B{S) 
such that 

1. T {A) and T (B) are disjoint whenever A, i? G B{S) are disjoint. 

2. T (^U°^i^j^ = UJL^T (Aj) whenever Ai,A2, • • • G B{S) are disjoint. 

Suppose /X is a measure on {S,B{S)) and let T be a full meta-time on S. Define /^i (• X T) as the 
mapping from B{S) into M given by /x A T) = ^ (T (A)) for any A G B{S). Then ^ (• A T) is 
a measure on {S, B{S)). We speak of /i (• A T) as the subordination of jj, by T. Similarly, if T is a 
meta-time on S we speak of /i (j4 A T) = /i (T (^4)) for A G Bb{S) as the subordination of /_i by T. 

Let us now recall an important result, which says that any measure, which is finite on compacts, 
can be represented as the image measure of the Lebesgue measure of a meta-time. 



Lemma 2. Bamdorff-Nielsen & Pedersen (2012 Lemma 3.1 ) Let T be a measure on B{S) satisfying 
T (A) < oo for all A G Bb{S) Then there exists a meta-time T such that, for all A G Bb{S), we have 

T{A) = Leb(T(A)). 

We speak of T as a meta-time associated to T. 
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Remark 9. T can be chosen so that T (A) = {T (x) : x G A} where T is the inverse of a measurable 
mapping U : 5 — )• 5 determined from T. Here integration under subordination satisfies 

/ f{s)f,{dsXT) = [ f{s)i2iT{ds))= [ f{U{u))fi{du). 
Js Js Js 

Suppose S is open and that T is a measure on S which is absolutely continuous with density r. If 
we can find a mapping T from 5 to sending Borel sets into Borel sets and such that the Jacobian of 
T exists and satisfies 

|T/xl = , 

then T given by T (A) = {T (x) : x G A}, is a natural choice of meta-time induced by T. In fact, by 
the change of variable formula, 

Leb{T{A)) = I lT(A)(y)dy = I Ia (x) |T/,|dx = ^ t (x) dx = T (A) , 

verifying that T is a meta-time associated to T. 

Remark 10. T is not uniquely determined by T and two different meta-times associated to T may 
yield different subordinations of one and the same measure /i. 

5.2.2 Extended subordination of Levy bases 

Let us now study the concept of subordination in the case where we deal with a Levy basis. Let T be 
a full random measure on S. Then, by Lemma |2j there exists a.s. a random meta-time T determined 
by T and with the property that Leb (T {A}) = T (A), for all A G Bb{S). There are two cases of 
particular interest to us. First, T is induced by a Levy basis L on 5 that is non-negative, dispersive 
and of finite variation. Second, T is induced by an absolutely continuous random measure T on 5 
with a non-negative density r satisfying T (^4) = X4 t (z) dz < 00 for all A G Bb{S). 

Definition 18. Let L be a Levy basis on S and let T, independent of L, be a full random measure on 
S. The extended subordination of L by T is the random measure L{- XT) defined by 

L{AXT) = L{T{A)) 

for all A£ Bb {S) and where T is a meta-time induced by T (in which case Leb (T (A)) = T (A)). 

Note that the Levy basis L may be n-dimensional. We shall occasionally write L AT for L (• XT), 
which is a random measure on {S, Bb{S)). 

In the case that L is a homogeneous Levy basis and since T is assumed to be independent of L, 
-L(- A T) is a (in general not homogeneous) Levy basis, whose conditional cumulant function satisfies 

C{C t L{A ^ T)\T} = T{A)C{CtL'} (38) 

for all A G Bb{S) and where L' is the Levy seed of L. On a distributional level one may, without 
attention to the full probabilistic definition of L(- X T) presented above, carry out many calculations 
purely from using the identity established in ( [38] ). 

Remarks 

The two key formulae L(- XT) = L{T (•)) and Leb (T {A)) = T (A) show that the concepts of ex- 
tended subordination and meta-time together generalise the classical subordination of Levy processes. 
Provided that L is homogeneous we have that 

C{C t L{A X T)\T} = Leb (T {A)) C{C \ L'}, 
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and hence 

C{C t L{A X T)} = logE {e^(^)C{aL'}| _ 
Hence we can deduce the following results. 

• The values of the subordination L{- XT) of L are infinitely divisible provided the values of T 
are infinitely divisible and L is homogeneous. 

• If L is homogeneous and if T is a homogeneous Levy basis then L(- A T) is a homogeneous 
Levy basis. 

• In general, T is not uniquely determined by T. Nevertheless, provided the Levy basis L is 
homogeneous the law of L{- X T) does not depend on the choice of meta-time T. 

Levy-Ito type representation of L (• A •) 

We have already reviewed the Levy-Ito representation for a dispersive Levy basis L, see Q. It follows 
directly that the subordination of L by the random measure T with associated meta-time T has a Levy- 
Ito type representation 

L (A A T) = a (T (A)) + ^/b(T{A})W (T {A)) 

+ [ yN {dy; T (A)) + [ y{N-v) {dy- T [A)) . 
J{\v\>l\ J\\v\<l\ 



Levy measure of L XT 



Suppose for simplicity that L A T is non-negative. According to Bamdorff-Nielsen (20101, the Levy 
measure D of L XT is related to the Levy measure i/ of T by 



i/(dx;s)= / P [L'Js) e dx) V {dy] s) 
Jr 



5.2.3 Extended subordination and volatility 

In the context of ambit stochastics one considers volatility fields a in space-time, typically specified 
by the squared field Tt (x) = (x). So far, we have used a stochastic volatility random field a in the 
integrand of the ambit field to introduce volatility modulation. 

A complementary method consists of introducing stochastic volatility by extended subordination. 
The volatility is incorporated in the modelling through a meta-time associated to the measure T on 5 
and given by 

T{A)= [ Tt (x) dxdt 
J A 

A natural choice of meta-time is T (A) = {T (x, t) : (x, t) G A}, where T is the mapping given by 

T(x,t) = (x,r+(x)) 
and where ^ 

(x) = -Ts (x) ds. 
Jo 

The above construction of a meta-time in a tempo-spatial model is very general. One can construct 
a variety of models for the random field rt(x), which lead to new model specifications. Essentially, 
this leads us back to the problem we tackled in the previous Subsection, where we discussed how such 
fields can be modelled. For instance, one could model Tt (x) by an Omstein-Uhlenbeck type volatility 



field or any other model discussed in Subsection 5.1 Clearly, the concrete choice of the model needs 



to be tailored to the particular application one has in mind. 
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5 TEMPO-SPATIAL STOCHASTIC VOLATILITY/INTERMITTENCY 



5.3 Probability mixing and Levy mixing 



Volatility modulation can also be obtained through probability mixing as well as Levy mixing. 

The main idea behind the concept of probability or distributional mixing is to construct new dis- 
tributions by randomising a parameter from a given parametric distribution. 

Example 19. Consider our very first base model ([T]l. Now suppose that the corresponding Levy basis 
is homogeneous and Gaussian, i.e. the corresponding Levy seed is given by V ~ N{fi + /Sa'^ , cr^) with 
/? G M and a"^ > 0. Now we use probability mixing and suppose that in fact a"^ is random. Hence, 
the conditional law of the Levy seed is given by L'\a ~ -/V(/i + cr^). Due to the scaling property 



of the Gaussian distribution, such a model can be represented as ( |16| ) and, hence, in this particular 
case probability mixing and stochastic volatility via a stationary stochastic integrand essentially have 
the same effect. Suppose that the conditional variance cr^ has a generalised inverse Gaussian (GIG) 
distribution rather than being just a constant, then L' follows in fact a generalised hyperbolic ( GH) 
distribution. Such a construction falls into the class of normal variance-mean mixtures. 

In this context it is important to note that probability /distributional mixing does not generally lead 
to infinitely divisible distributions, see e.g. Steutel & van Ham ( |2004[ Chapter VI). Hence Bamdorff- 
Nielsen, Perez-Abreu & Thorbj0msen (2012 1 propose to work with Levy mixing instead of probabil- 
ity/distributional mixing. Levy mixing is a method which (under mild conditions) leads to classes of 
infinitely divisible distributions again. Let us review the main idea behind that concept. 

Let L denote a factorisable Levy basis on with CQ (a, b, v{dx),c). Suppose that the Levy 
measure v{dx) depends on a possibly multivariate parameter G 0, say, where denotes the pa- 
rameter space. In that case, we write v{dx; 6). Then, the generalised Levy measure of L is given by 
v{dx; 6)c{dz). Now let 7 denote a measure on and define 



where we assume that 



h(dx,dz) = / u{dx;9)j{d0)c{dz) 
Je 



/(lA.= ),X<i.,<iz)<oo. 

Jr 



(39) 



Then there exists a Levy basis L which has h as its generalised Levy measure. We call the Levy basis 
L the Levy basis obtained by Levy-mixing L with the measure 7. 
Let us study a concrete example of Levy mixing in the following. 

Example 20. Suppose L is a homogeneous Levy basis with CQ given by (0,0, i'{dx), Leb) with 
v{dx) = 98i{dx) for 9 > 0. I.e. the corresponding Levy seed is given by L' ~ Poi{6). Now we 
do a Levy-mixing of the intensity parameter 6. Let 



n{dx, dz) 



ej{d9)5i{dx)dz, 



for a measure 7 satisfying condition ( 39 1. Let L be the Levy basis with CQ (0, 0, /q 6j{dv)5i (dx) , Leb). 
In that case, the base model ([T]l would be transferred into a model of the form 

I I h{:x.,t;^,s)L{d^,ds,dv). 

Je JAtiyi) 

Example 21. Let us consider the example of a (sup)OU process, see Barndorff-Nielsen ( |2000 1, 



Bamdorff-Nielsen & Stelzey\^201I\ and \BarndoTff-Nielsen, Perez-Abreu & Thorbj0msen\^2OI2\ . Let 

-'^'-'Ul{s), 



L denote a subordinator with Levy measure ul (and without drift) and consider an OU process 

ft 

Yt 



> 0. 
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6 CONCLUSION AND OUTLOOK 



A straightforward computation leads to the following expression for its cumulant function (for C, G Mj.' 

[e'^"" -l^v{dx;e), where u{dx,e) = j VL{e'^''dx)du 

is a mixture of ul with the Lebesgue measure. A Levy mixing can be carried out with respect to the 
parameter 9, i.e. 

PCX) 

i){dx) = / u{dx,9)j{de), 
Jo 

where 7 is a measure on [0, 00) satisfying xD(dx) < 00. Then u is a Levy measure again. Now, 
let L be the Levy basis with extended Levy measure 

VL{dx)du'y{d9), 

and define the supOU process Yt w.rt. L by 

f'OO ft 

Yt= / e-^(*-^)L((is,(iA). 

JO J -oa 

Then the cumulant function ofYt is given by 

/* 00 /*i /*oo /* 00 

C{C \Yt}= / C{Ce-^(*-^) \ Li}dsj{d9) = / / (e^^^"'" - l) UL{dx)du-f{d9) 
Jo J-00 Jo Jo ^ ' 

= [e'^"" - 1) D{dx). 
Hence, we have seen that a supOU process can be obtained from an OU process through Levy-mixing. 
5.4 Outlook on volatility estimation 

Once a model is formulated and data are available the question of assessment of volatility arises and 
while we do not discuss this in any detail, tools for this are available for some special classes of ambit 
processes, see |BamdorfF-Nielsen et al.| ( [2009[ ), |BamdorfF-Nielsen, Corcuera & Podolskij| ( |2011[[20T2] ) 



and also Bamdorff-Nielsen & Graversen (2011 1. Extending these results to general ambit fields and 



processes is an interesting direction for future research. 

6 Conclusion and outlook 

In this paper, we have given an overview of some of the main findings in ambit stochastics up to 
date, including a suitable stochastic integration theory, and have established new results on general 
properties of ambit field. The new results include sufficient conditions which ensure the smoothness 
of ambit fields. Also, we have formulated sufficient conditions which guarantee that an ambit field is 
a semimartingale in the temporal domain. Moreover, the concept of tempo-spatial stochastic volatil- 
ity/intermittency within ambit fields has been further developed. Here our focus has been on four 
methods for volatility modulation: Stochastic scaling, stochastic time change and extended subordi- 
nation of random measures, and probability and Levy mixing of the volatility/intensity parameter. 

Future research will focus on applications of the general classes of models developed in this paper 
in various fields, including empirical research on turbulence modelling as well as modelling e.g. the 
term structure of interest rates in finance by ambit fields. In this context, it will be important to 
establish a suitable estimation theory for general ambit fields as well as inference techniques. 
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